Ab initio study of 2p-core level x-ray photo-emission spectra in ferromagnetic 

transition metals 
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We study the 2p-core level x-ray photo-emission spectra in ferromagnetic transition metals, Fe, Co, 
and Ni using a recently developed ab initio method. The excited final states are set up by distributing 
electrons on the one-electron states calculated under the fully screened potential in the presence of 
the core hole. We evaluate the overlap between these excited states and the ground state by using 
one-electron wave functions, and obtain the spectral curves as a function of binding energy. The 
calculated spectra reproduce well the observed spectra displaying interesting dependence on the 
element and on the spin of the removed core electron. The origin of the spectral shapes is elucidated 
in terms of the one-electron states screening the core hole. The magnetic splitting of the threshold 
energy is also estimated by using the coherent potential approximation within the fully screened 
potential approximation. It decreases more rapidly than the local spin moment with moving from 
Fe to Ni. It is estimated to be almost zero for Ni despite the definite local moment about 0.6 fj,s, in 
agreement with the experiment. 

PACS numbers: 79.60.-i 71.15.Qe 71.20.Be 
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I. INTRODUCTION 

Core-level x-ray photo-emission spectroscopy (XPS) is 
one of the powerful tools for studying the electronic struc- 
ture in solids through the response of electrons to the 
photocreated core hole. It is well known that the re- 
sponse function in metals displays singular behavior near 
the Fermi edgei~— Core-level XPS spectra as a func- 
tion of binding energy display asymmetric shape near 
the threshold. 1 * The spectra sometimes have extra struc- 
tures in the high binding energy region. A notable ex- 
ample in metals is a satellite peak on the 2p XPS in 
Ni metal, which appears around 6 eV higher than the 
threshold^ Feldkamp and Davis& analyzed the Ni 2p XPS 
spectra by evaluating the overlap determinants between 
the ground and excited states, using a numerical method 
on the Hubbard-like model. They clarified the origin of 
satellite as a combined effect of the core-hole screening 
and the interaction between electrons, and estimated the 
strength of the core-hole potential consistent with the 
binding energy of the satellite intensity and the asymme- 
try parameter. 

In addition to the above features, core-level XPS spec- 
tra present several intriguing behaviors in ferromagnetic 
transition metals Fe, Co, and Ni. On the spin resolved 
3s spectra, they show characteristic satellite intensities 
in the majority spin channel for all Fe, Co, and Ni, while 
they show almost single peak structures in the minority 
spin channel. In Fe, the spectra exhibit particularly a 
large satellite structure only in the majority spin chan- 
nel. In contrast to the 3s spectra, the 2p spectra do 
not have a clear satellite peak in Fe, while they show 
the notable 6-eV satellite in Ni. As regards the peak 
around the threshold in 2s, 2p, and 3s spectra, its posi- 
tions noticeably depend on the spin channel in Fe, while 
such dependence has not been observed in Niii& Accord- 
ing to the recent hard x-ray photoemission spectroscopy 



(HAXPES) experiment, 9 such magnetic splitting is esti- 
mated about 0.9 eV for the 2p 3 / 2 core in ferromagnetic 
Fe. Similar magnetic splittings have also been observed 
on the 2p 3 /2 spectra in several half-metallic ferromag- 
netic Heusler alloys J£r— Such splittings are usually con- 
sidered to be related with the local spin moment at the 
photo-excited site. As regards the satellite structure, it 
has been clearly observed in the 3s XPS in several Fe 
compounds. 7 We have clarified the origin of the satellite 
intensity on the 3s spectra and that the satellite inten- 
sity is not a direct reflection of the local spin moment by 
calculating the XPS spectra on the ab initio leveh 15 i 16 

For the strongly correlated localized electron systems 
such as the 3d transition metal oxides and the /-electron 
systems, the theoretical analysis of the XPS spectra have 
been carried out mainly on the basis of atomic, cluster, 
or impurity Anderson models^ Although the impurity 
Anderson model has been applied to analyze the core- 
level spectra in some itinerant metallic systems; 12 ! 18 it 
is not suitable to analyze the spectra in the highly itin- 
erant metallic systems. Another approach to investigate 
the core-level spectra is based on the independent one- 
electron theory exploiting a meanfield approximation. 
Mahan,— and Barth and Grossmann22 have calculated 
the overall line shapes of the x-ray spectra by evaluating 
one-electron wave functions under the final state poten- 
tial, and have reproduced well the experimental emission 
and absorption spectra in metals within the meanfield 
theory. 

Recently we have presented the XPS spectra based on 
the ab initio band structure calculation! 15 i 16 We have cal- 
culated the final state potential self-consistently by car- 
rying out the band structure calculation on the system 
holding a core- hole at a photo-excited site. Not only the 
core-hole potential but also the relaxation of the core 
states as well as the screening electron distribution are 
automatically determined through the calculation. Dis- 



tributing electrons on the energy levels thus evaluated, we 
have constructed various final states. We have evaluated 
the overlaps between those final states and the ground 
states by using the one-electron wave functions, and fi- 
nally obtained the 3s core-level XPS spectra in the ferro- 
magnetic transition metals Fe, Co, and Ni. The element 
and photo-electron-spin dependence of the spectral line 
shape have been reproduced in good agreement with ex- 
periments within the ab initio manner ^ 5 ' 16 This method 
may be regarded as an extension of the Feldkamp-Davis 
method* to an ab initio level. The origin of the spectral 
shape as a function of the binding energy has been elu- 
cidated in terms of the one-electron states screening the 
core-hole as follows. 

In the fully screened state where an up-spin 3s elec- 
tron is removed, the down-spin 3d states are strongly at- 
tracted in the core-hole site, forming quasibound states 
near the bottom of the 3d band (hereafter majority spin 
and minority spin are called as up-spin and down-spin, 
respectively). The final states that the quasibound states 
are unoccupied contribute to the satellite or shoulder in- 
tensities. In the state where a down-spin 3s electron is 
removed, the up-spin 3d electrons are strongly attracted 
to the core hole, forming quasibound states near the bot- 
tom of the 3d band, but the down-spin 3d electrons are 
not attracted strongly enough to form quasibound states. 
Because the up-spin 3d states are almost fully occupied 
in the ground state, the quasibound states, which appear 
only in the up-spin 3d states, could not become unoccu- 
pied, leading to a single peak structure of the spectra. 

In this paper, we apply our method to calculate the 2p 
core-level XPS spectra in ferromagnetic Fe, Co, and Ni, 
and elucidate the underlying mechanism. The calcula- 
tion could not distinguish between the 2p 3 / 2 core and the 
2pi/2 core, because the spin-orbit interaction is not taken 
into account. We focus on the difference of the spectra 
among the three ferromagnetic transition metals and on 
the difference between the 2p and 3s spectra. The main 
features of the 2p spectra are consistently reproduced in 
the same manner as those of 3s spectra. We could under- 
stand the origin of the spectral shape in the high binding 
energy region by considering an electron-hole pair exci- 
tation from the quasibound states with down-spin to the 
unoccupied states in the fully screened state. In Ni, the 
fully screened one-electron states are similar both in the 
2p and 3s electron removal states, while they are signifi- 
cantly different in Fe. Such different screening behaviors 
of the 3d states cause the difference in the spectral shapes 
between Fe and Ni. It is found that the band filling in the 
down-spin state plays important roles to give rise to the 
difference between the Fe 2p and Ni 2p spectra. We also 
evaluate the magnetic splitting of the threshold energy, 
exploiting the coherent potential approximation (CPA) 
within the fully screened potential approximation. This 
may be observed as the magnetic splitting of the peak 
position around the threshold in the spin-resolved XPS 
spectra. We obtain the magnetic splitting of Fe 2p 3 / 2 as 
0.9 eV in agreement with the recent experiment, and that 



of Ni 2p 3 / 2 as nearly zero despite the finite local magnetic 
moment about 0.6 ^b- 

The present paper is organized as follows. In Sec. II, 
we briefly describe the calculation procedure of the XPS 
spectra. In Sec. Ill, we present the 3d band calculated in 
the presence of core hole and the XPS spectra in compar- 
ison with the experiments. The last section is devoted to 
the concluding remarks. 



II. CALCULATION METHOD 
A. Calculation of spectral intensity 

We briefly summarize the calculation procedure for the 
core-level photoemission spectra. For details, we refer the 
readers to Ref. [lg| . The many-body wave functions of the 
ground and final states are assumed to be given as single 
Slater determinants. 

For the ground state, we carry out the full potential lin- 
ear augmented plane wave (FLAPW) band structure cal- 
culation based on the local density approximation (LDA), 
and obtain the one-electron wave functions 0j's. We con- 
struct the Slater determinant by putting electrons from 
the lowest energy state up to the Fermi level. For the 
final states, we carry out the same type of band calcu- 
lation under the fully screened potential in a periodic 
array of supercells with one core hole per cell. The core- 
states are treated as localized states within a muffin-tin 
sphere, so that we could specify the core-hole site. In 
reality, only one core-hole should exist in crystal through 
the XPS event. Therefore, the larger the unit cell size is, 
the better results would be expected to come out. We 
use the 3x3x3 bcc supercell for Fe and fee supercell 
for Co and Ni, where the core-hole sites form a bcc super 
lattice and an fee super lattice, respectively. In the self- 
consistent procedure, we keep a hole in a specified core 
level at a core-hole site, and put an extra electron in each 
super cell to guarantee the charge neutrality. The local 
charge neutrality would be satisfied, known as the Friedel 
sum rule in the impurity problem, 21 We thus obtain the 
one-electron wave function ipi with energy eigenvalue q, 
which takes account of not only the effect of core-hole 
potential but also that of electron-electron interactions 
within the limit of the LDA. The final state |/ (0)) with 
the lowest energy is constructed by putting electrons from 
the lowest energy state up to the Fermi level at each k 
point, as was done in the ground state. The other final 
states are obtained by creating electron-hole (e-h) pairs 
from this final state \f (0)}. We designate the state hav- 
ing the I e-h pairs on the state \f (0)) as \f (£, m)), where 
the index m distinguishes the different electron config- 
uration. The overlap integral a\ between the wave 
function i/y^m) for the z'th occupied valence state in the 
final state \f (£, m)) and the wave function <\>j for the j'th 
occupied valence state in the ground state \g) is given by 

1 = f ipin m^jdr 3 , where the volume integral is 



as a,- 



i(e,m)^3 l 



taken over a super cell. 

Neglecting the interaction between the escaping photo- 
electron and the other electrons in matter, we consider 
the XPS process that a core-electron is excited to a high 
energy state with energy e by absorbing an x-ray photon 
with energy lu. Exploiting Fermi's golden rule, we obtain 
the expression of the spectral intensity as a function of 
the biding energy u> — e as 
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where Eq and E g represent the total energy of the fi- 
nal state |/(0)) and the ground state \g) , A is an en- 
ergy independent constant including the contribution of 
overlaps between the wave functions of the core-state in 
the final and ground states. iV c is the number of va- 
lence electrons in the ground state. The overlaps between 
the valence states and the core states are eliminated be- 
cause they could be almost orthogonal. A£^ m ) is the 
excitation energy defined by AEr^ tm ) = Eii >m -\ — Eq — 
X\ n')( e « — e n '), where E(t,m) represents the total en- 
ergy of the final state \f (£, m)) and e n 's are the Kohn- 
Sham eigenvalues. The energy difference e„ — e n i stands 
for the energy of an e-h pair of an electron at level n 
and a hole at level n' and the summation are taken over 
all e-h pairs in the final state \f(£,m)). Although the 
Kohn-Sham eigenvalues may not be proper quasi parti- 
cle energies, they practically give a good approximation 
to quasiparticle energies, except for the fundamental en- 
ergy gapj 22 ' 23 In the following calculation, we replace the 
5 function by the Lorentzian function with the full width 
at half maximum (FWHM) 2T = 1 eV in order to take 
into account the lifetime broadening of the core level. 
Eq — Eq is treated as an adjustable parameter so that 
the threshold of XPS spectra coincides with the experi- 
ments. In order to suppress the error caused by the ficti- 
tious periodicity of the core-hole site, we pick up only the 
r point as the sample states for calculating XPS spectra. 
For Ni, we pick up the X point (i, 0, 0) as a sample point 
in addition to the T point, because the 3d band states at 
the r point are fully occupied by both up- and down-spin 
electrons on the system of the 3 x 3 x 3 fee super-cell. We 
take account of the final states including 0, 1, 2, and 3 
electron- hole pairs on the final state \f (0)) and restrict 
the final states \f (£, m)) so that the excitation energy 



E, 



Eq is less than 10 eV. We need to calculate the 



more than 10 determinants of the matrices of the size of 
about 150 x 150 for Fe even after the above simplification. 
Before closing this subsection, we briefly mention the 
limitations in this calculation. First, we assume that 
the 2p core hole is spherically distributed with neglect- 
ing the dependence on the magnetic quantum number 
of the core hole. Second, we take no account of the 
spin-orbit interaction (SOI) in the band structure cal- 



culation. Third, due to the finiteness of the cell size, 
the final state |/(0)) (no e-h pair) has a finite overlap 
with the ground state \g), resulting in a finite intensity 
at the threshold. In principle, such overlap should con- 
verge to zero with N e — > oo, according to the Anderson 
orthogonality theorem^ energy levels become continuous 
near the Fermi level and thereby infinite numbers of e-h 
pairs could be created with infinitesimal excitation en- 
ergies, leading to the so-called Fermi edge singularity in 
the XPS spectra. The finite contribution obtained above 
arises from the discreteness of energy levels and could be 
interpreted as the integrated intensity of singular spectra 
near the threshold, in consistent with the model calcula- 
tions for other systems w^ 



B. Binding energy difference at the threshold 



The threshold energies wi 11 for the up-spin core- 
electron removal excitation and tuf 1 for the down-spin 
one are generally different in the ferromagnetic systems. 
It could be naively considered that the difference Aw th = 
wi — u;j h is linked to the energy difference between the 
up- and down-spin core levels in the ground state and 
is a good indicator for the local spin moment. However, 
because of the considerable core-hole screening in met- 
als, the final state effects should be taken into considera- 
tion. Within the fully screened potential approximation, 
we may obtain a better estimate of Aw th without using 
supercell but by using the CPA in the low concentra- 
tion limit, as was done in the calculation of the core-level 
chemical shift in metallic random alloysj 25 ' 26 Because the 
core-hole screening is almost completed within the core- 
hole site for metals, the CPA, which is a one-site approx- 
imation, may not cause large error. 

The threshold energy may be written as 



, ,th _ j? 
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where E g represents the total energy of the ground 
state, -Ej-m does that of the fully screened state with 
an up(down)-spin core hole and N e + 1 band electrons, 
and £f is the Fermi energy in the fully screened state. 
Therefore, the difference may be given by 
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The energy difference £Wn — E g corresponds to the en- 
ergy of putting one core-ionized impurity at a normal site. 
This might be given as the generalized thermodynamic 
chemical potential, 
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where E-^n ) represents the total energy per unit cell for 
the random alloy system which consist of core-ionized 
and normal atoms and p the density of the cell including a 



core-hole site. Hence, the difference Aw th may be written 
as 



Au 



th 



dp 
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We evaluate these values with the help of the KKR 
Green's function band structure calculation combined 
with CPA.— 12£ We carry out the calculation at the con- 
centration ranged from p = to 0.05. Aw th is calculated 
by an extrapolation to zero concentration p. 



III. RESULTS AND DISCUSSIONS 

A. Ground and fully screened states 

1. d-DOS 

We carry out the band structure calculation based 
on the LDA. For the 2p electron removal, we make the 
FLAPW super cell calculation. The magnetic moments 
in the ground states are obtained 2.1, 1.6, and 0.6 pb in- 
side a muffin-tin sphere (r m = 2.0 Bohr) for Fe, Co, and 
Ni, respectively. The obtained electronic structures in 
the ground states are consistent with the past band struc- 
ture calculations^ We have also checked the accuracy of 
the calculation by carrying out the KKR band structure 
calculation combined with the CPA j 27 ' 28 We treat the 
system as an random alloy system consisting of the nor- 
mal and core-ionized atoms, and take the limit of zero 
concentration of the core-ionized atom. The FLAPW 
supercell calculation and the KKR-CPA calculation give 
essentially the same results. 

The one-electron states in the fully screened state are 
modified from those in the ground state by the core-hole 
potential. Such changes in one-electron states may be ex- 
plained through the change in the d-symmetric density 
state (d-DOS) at the core-hole site. Figure Q] shows the 
d-DOS at the core-hole site. The d-DOS's at the core- 
hole site in the fully screened states are greatly modified 
from those in the ground state. The d-DOS at sites with 
no core hole, which are not shown here, are almost un- 
changed from those in the ground state. This indicates 
that the core-hole potential is almost completely screened 
inside the core-hole site, being consistent with the small 
screening length ~Bohr radius for metals. The change of 
the d-DOS's by the core-hole potential depends largely 
on whether the up- or down-spin electron is removed from 
the 2p G core state. We refer to the hole created in the 
2p states by removing an up(down)-spin electron as the 
up(down)-spin hole. 

In the fully screened state with an up-spin 2p hole, 
the weight of the down-spin d-DOS moves toward the 
higher binding energy region if compared with the up- 
spin d-DOS. This is due to the repulsive (attractive) ex- 
change interaction working between the up-spin hole and 
the up(down)-spin 3d electrons, in addition to the attrac- 
tive interaction between the core hole and 3d electrons. 



In Ni, the weight of the d-DOS of the down-spin state 
concentrate near the bottom of the 3d band, indicative of 
the quasibound states. The weight of the up-spin states 
also shifts toward to the higher binding energy region. 
These behaviors are very similar to those for the 3s elec- 
tron removal. In Co, the weight of the d-DOS's of the 
down-spin 3d states shifts toward the higher binding en- 
ergy region. Opposite to the case for Ni, that weight 
of the up-spin 3d states shifts toward the Fermi level. 
The behaviors are similar to those for the 3s electron re- 
moval but with rather moderate shape. In Fe, the weight 
of the d-DOS's of the down-spin 3d state moves toward 
the higher binding energy region, but not as much as for 
Ni. The weight of the up-spin 3d states is almost un- 
changed. The behavior in Fe contrasts to that for the 3s 
electron removal, where the weight of the down-spin d- 
DOS highly concentrates near the bottom of the 3d band 
and that for the up-spin d-DOS significantly shift toward 
the Fermi level. 

In the fully screened state with a down-spin 2p hole, 
the weight of the up-spin d-DOS's tends to move toward 
the higher binding energy region probably due to the 
exchange interaction between the core hole and the 3d 
electrons. In Fe, the weight of the up d-DOS is mod- 
erately shifted toward the higher binding energy region 
compared to the case of the down-spin 3s electron re- 
moval. The weight of the down-spin d-DOS also slightly 
moves toward the higher binding energy region. In Ni, 
the shift of the up-spin d-DOS is large, indicative of the 
quasibound states. The weight of the down-spin d-DOS 
also concentrates near the bottom of the 3d band. This 
behavior in Ni is very similar to that for the down-spin 
3s electron removal. The behavior of the d-DOS's in Co 
is intermediate between Fe and Ni. To sum up, the 3d 
states in Fe for the 2p electron removal are quite different 
from those for the 3s electron removal, while those in Ni 
are quite similar in both removal states. 



2. Screening electron number 

The change of the d-DOS's is related to the screen- 
ing electron number at the core-hole site. Table U shows 
the screening electron numbers Aridf, An^i, and And-|- + 
An dl . Here An dt(i) is given by An dt{i) 



where 



dt(l) 



MU) n dt(4)' 

is the up(down)-spin electron number inside 
the muffin-tiii sphere in the d-symmetric states in the 
core-hole site in the fully screened state, and n fL ,-, is the 
corresponding quantity in the ground state. Total screen- 
ing electron numbers Arid-]- + An^ inside the muffin-tin 
sphere is found to be roughly unity. On the sites with 
no core hole, the electron numbers are almost unchanged 
from that in the ground state, indicating that the core- 
hole screening is almost completed inside the core-hole 
site. 

First we discuss the case that an up-spin 2p electron 
is removed. In the fully screened state of Co, the core 
hole attracts so strongly the down-spin 3d electrons that 
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Figure 1: Density of state projected onto the state with d symmetry inside the muffin-tin sphere (d-DOS) at the core-hole site 
in Fe, Co, and Ni. Top panels show the d-DOS's in the ground states. Middle (bottom) panels show those in the fully screened 
states with an up (down)-spin 2p core-hole. The d-DOS's in the fully screened states with a 3s hole are also shown as inset. 
The energy of the Fermi level is zero. 



they overscreen the core hole, and, as a countereffect, a 
few up-spin 3d electrons at the core-hole site are pushed 
away from the muffin-tin sphere due to the Coulomb re- 
pulsion with the excess down-spin electrons. In the fully 
screened state of Ni, although the strong attraction by 
the core hole is expected to work on the down-spin elec- 
trons, the screening down-spin electron number is less 
than unity, because the d states at the core-hole site are 
already almost filled in the ground state. The up-spin 
3d electrons participate a little to screen the core-hole. 
Note that those screening electron numbers are similar 
to those for the up-spin 3s electron removal state both 
in Co and Ni. In Fe, the screening of the core hole is 
almost completed by the down-spin 3d electrons with- 
out overscreening. This is quite different from the fully 
screened state for the up-spin 3s electron removal, where 
the down-spin 3d electrons extremely overscreen the core 
hole, and that overscreening is compensated by pushing 
away the up-spin 3d electrons from the core hole site. 
Now we discuss the case that a down-spin 2p electron is 
removed. In the fully screened states of Ni, although the 
up-spin 3d electrons are to be strongly attracted by the 
core hole, the screening electron number is not large, be- 
cause the up-spin 3d state is almost filled in the ground 
state. Therefore, the down-spin 3d electrons contribute 



Table I: Screening electron number with the d symmetry in- 
side the muffin-tin sphere at the 2p core hole site. The radii 
of the muffin-tin spheres are 2.0 Bohr. 
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0.19 
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0.67 


1.14 


Co 


Up 


-0.14 
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0.22 
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1.08 


Ni 


Up 


0.18 


0.89 


1.07 




Down 


0.27 


0.79 


1.06 



in large amount, leading that the total screening num- 
bers become nearly unity. In the fully screened states of 
Fe and Co, although the up-spin 3d electrons are more 
attracted by the core hole than the down-spin electrons, 
the screening electron number of the up-spin electron is 
smaller than the down-spin electrons for a similar reason 
to Ni. 



3. Difference between the 2p core electron removal and the 
3s-core electron removal 



Now we focus on the difference in the screening behav- 
ior between with the 2p core hole and the 3s core hole. 
For both cases, the down(up)-spin 3d electrons are, gen- 
erally speaking, attracted to the core hole more strongly 
than the up(down)-spin 3d electrons in the fully screened 
state both with an up(down)-spin 2p or 3s core hole, be- 
cause the exchange interaction is working between the 
core hole and the 3d electrons. In addition to this ten- 
dency, the screening in Fe seems moderate in the fully 
screened state with a 2p core hole, compared to the pres- 
ence of the overscreening with a 3s core hole. On the 
other hand, in Ni, the screening with a 2p core hole is 
comparably larger than that with a 3s core hole. How 
does such a different screening behavior occur in Fe and 
Ni? The exchange interaction J^p-zd between the 2p and 
3d electrons, and Jz s -3d between the 3s and 3d electron 
may be estimated by the atomic Hartree-Fock calcula- 
tion, which are J2p-3d ~ 1 eV and J^s-zd ~ 2 eV, for 
all the neutral Fe, Co, and Ni atoms, respectively. The 
small 2p-3d exchange interaction in Fe could not solve 
the above question. 

To eliminate the effect of the 2p-3d (3s-3d) exchange 
interaction, we carry out the band calculation under the 
condition that a half up-spin and a half down-spin elec- 
tron is removed from the 2p (3s) core states. Figure 
[5] shows the d-DOS's at the core hole site in the fully 
screened state. The d-DOS's for the state with a 2p core 
hole are similar to those with a 3s core hole. Note that 
the atomic Hartree-Fock calculation gives quite different 
values of the Slater integrals F°'s for 2p-3d and 3s-3d in- 
teractions, that is, F° (2p, 3d)-F° (3s, 3d) ~ 7eV, which 
in fact contradicts the results shown in Fig. [5J In Ni, 
the d-DOS's concentrate near the bottom of the up- and 
down-spin 3d band, forming quasibound states. In Fe, 
the situation is quite different; the weights of the d-DOS 
shift only moderately toward the higher binding energy 
region, almost keeping the structures of the ground state. 
In Co, the behavior of the d-DOS is intermediate between 
Fe and Ni. This difference in Fe, Co, and Ni may be 
related to the different electron occupation in the band 
in the ground state. In Ni, because the 3d band is al- 
most fully occupied, it is hard to make enough room for 
accommodating the screening electron at the core hole 
site by mixing the one-electron states within the 3d band 
states. Therefore, the mixing beyond the 3d band states 
has to take place to complete core hole screening inside 
the core hole site. These processes might help to form 
quasibound states. In Fe, because the 3d band in the 
ground state has enough room to accommodate excess 
screening electron, the mixing of the one-electron states 
within the 3d band states might be sufficient in order to 
complete core hole screening, resulting in only the slight 
shift of the weight of the d-DOS's. When the 2p-3d (3s- 
3d) exchange potential is turned on, the screened states 
get modified further. In Ni, both the 2p-3d and 3s-3d ex- 



change potentials simply shift the quasibound states to 
the higher or lower binding energy regions. In Fe, the 3s- 
3d exchange potential considerably modify the 3d state 
by pulling down the down-spin 3d states at around 0.5 eV 
above the Fermi level (the top panel in Fig. [5]) under the 
Fermi level, and by making the quasibound states formed 
near the bottom of the down-spin 3d band. The 2p-3d 
exchange potential, on the other hand, hardly affects the 
3d band states, that is, the down-spin 3d states at around 
0.5 eV still stay above the Fermi level probably due to 
the smallness of the 2p-3d exchange potential. 



B. Photoemission Spectra 

Using the one-electron wave functions in the ground 
state and fully screened states, the XPS intensities are 
calculated as a function of the binding energy from eq. 
(UJ). Figures [5H5] show the calculated spectra for bcc Fe, 
fee Co, and fee Ni, in comparison with the experimen- 
tal 2p 3 /2 spectra. Although the edge singularity is not 
reproduced, the overlaps (/ (0) \g) and (/ (l,m) \g)'s ap- 
pear to give the reasonable intensities around threshold. 
The total spectral curves, which depend on the elements, 
reproduce well the overall structures observed in the ex- 
periments. Note that the final states \f (I, m))'s holding 
a down-spin e-h pair mainly contribute to the intensity. 
The final state |/(l,m))'s created by putting an up e- 
h pair on the final state \f (0)) hardly contribute to the 
intensity, because the up-spin 3d band states are almost 
filled in the ground state, and thereby the overlap deter- 
minant become quite small. 

The Fe spectra show a single-peak structure with a 
weak shoulder intensity for both the up- and down-spin 
2p electron removals, contrasting with the strong satellite 
intensity for the up-spin 3s electron removal. Because the 
one-electron wave functions in the fully screened state are 
not so strongly modified from those in the ground state, 
the overlap (/ (0) \g) is almost unity. Consequently, the 
other overlaps (/ (£, m) |g)'s are almost zero, because the 
one-electron wave functions for the unoccupied states in 
the final state \f (0)) are nearly orthogonal to those for 
occupied states in the ground state \g). Note that the ob- 
served 0.9 eV splitting of the peak position around the Fe 
2p 3 / 2 threshold is not reproduced by the e-h excitations 
on the final state |/(0)). We expect that the splitting 
originates from the different threshold energy for the up- 
and down-spin 2p electron removals. 

In Co, the intensities at the higher binding energy re- 
gion are larger than the Fe 2p spectra, forming a rather 
strong shoulderlike structure for the up-spin 2p electron 
removal. Because the one-electron wave functions for the 
down-spin 3d band states in the fully screened state are 
noticeably modified from those in the ground state, the 
one-electron wave functions of the unoccupied levels in 
the final state \f (0)) have some amplitudes of those of 
the occupied levels in the ground state. The overlaps 
(/ (£, m) |<jf}'s accordingly become finite, giving rise to the 
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Figure 2: d-DOS's when removing a half up- and a half down-spin electron from the 2p and 3s core. The energy of the Fermi 
level is zero. 



shoulder. For the down-spin 2p electron removal, the 
spectra show a single-peak structure, because the down- 
spin 3d states are not strongly modified by the down-spin 
2p hole. 

The Ni spectra exhibit satellite structure for the up- 
spin 2p electron removal. The e-h pair excitations that an 
electron on the quasibound state is excited to the unoccu- 
pied states give rise to the satellite intensity. The satellite 
peak position is around 4 eV, which is 2 eV smaller than 
the so-called 6-eV satellite observed in experiment. This 
discrepancy might be caused by the LDA or the fully 
screened potential approximation. Because both the up- 
and down-spin 3d states are fully occupied in the fully 
screened state, the 3d-3d Coulomb interaction may not 
be relevant to this discrepancy within the fully screened 
potential approximation. Braicovich and van der Laan 32 
estimated the screening time constant r s = 1.5 fs in Ni, 
which is definitely longer than the values 0.18 fs in Fe 
and 0.43 fs in Co. Because the constant t s is comparable 
to the core hole lifetime r c ~ 1 fs, the use of the fully 
screened potential may not be appropriate to describe 
the satellite intensities in Ni. The spectra for the down- 
spin electron removal also show relatively large shoulder 
intensities, which would come from the excitations in the 
down-spin 3d states, because the unoccupied levels are 
available in the down-spin 3d states in the fully screened 
3d states. 



C. Energy difference of the threshold 

In order to estimate the magnetic splitting of the 
threshold energy, we calculate the energy difference of 
the threshold in the up- and down-spin core-electron re- 
movals by using the KKR-CPA method, as discussed in 
Scc lHBl Figures [5] and [7] show the energy difference of 
the threshold Aw th defined by Eq. (JSJ and the energy 
difference of the core levels Ae c — e? — el in the ground 



states. While Ae c 's are nearly proportional to the lo- 
cal spin moment, Aw th 's decrease more rapidly than the 
change of the local moments with moving from Fe to Ni. 
In Fe, we get the local moment about 2.2 /xb and Aw th 's 
about 1.2, 0.9, and 1.5 eV for the 2s, 2p, and 3s exci- 
tations, respectively. In Ni, on the other hand, Aw th 's 
are almost zero, despite the definite local spin moment 
about 0.6 /ib- These values seem to be consistent with 
the experimental observations; the splitting is estimated 
as > 1 eV in the Fe 3s spin resolved spectra^— and 0.9 
eV in the Fe 2p 3 / 2 HAXPES spectra. 9 On the other hand, 
such splittings have not been observed in the Ni 2p and 
3s spectra. 

It is obvious that the core hole screening plays crucial 
roles in determining the magnitude of the splitting of the 
threshold energy. The fact that the magnetic splitting 
in Ni is smaller than that in Fe may be understood as 
follows. All the 3d states in the fully screened state in 
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Figure 3: Fe 2p XPS spectra. Gray thick curve represents 
the total intensity. The curves denoted / , up /' , dn p > , 
/ , and _p 3 ' are the contribution from the final states in- 
cluding zero, one up-spin e-h pair, one down-spin e-h pair, 
two e-h pairs, and three e-h pairs on the final state |/(0)), 
respectively. Experimental data are taken from Ref. 0, which 
are not spin resolved. 



Figure 4: Co 2p XPS spectra.Gray thick curve represents the 
total intensity. The curves denoted / (0) , up / (1) , dn / (1) , / (2) , 
and /' 3 ' are the contribution from the final states including 
zero, one up-spin e-h pair, one down-spin e-h pair, two e-h 
pairs, and three e-h pairs on the final state |/(0)}, respec- 
tively. Spin-resolved experimental data are taken from Ref. 

M 



Ni are are pulled down below the Fermi level forming the 
quasibound states resulting in a suppression of the local 
spin moment at the core hole site solely by the Coulomb 
interaction between the core hole and the 3d electrons 
(Fig. [2]). The exchange potential could not give rise to a 
large difference in the screening electron distribution be- 
tween for the up-spin and down-spin electron removals. 
The 3d states in Fe, on the other hand, are strongly af- 
fected by the exchange potential particularly for the 3s 
electron removal, with the considerable change of the d- 
DOS's and with the screening electron numbers. Even 
for the 2p electron removal, the screening electron num- 
bers depend considerably on the spin channel. Thus we 
may guess that the weak effect of the exchange potential 
on the 3d states in Ni may be the origin of the magnetic 
splitting of the threshold energy. Note that the zero mag- 
netic splitting does not necessarily mean that the XPS 
spectra for the up-spin channel are the same as those for 
the down-spin channel. 



IV. CONCLUDING REMARKS 

We have applied an ab initio method to calculate the 
2p core-level XPS spectra as a function of binding energy 
in ferromagnetic metals Fe, Co, and Ni. The calculated 



spectra have been compared to the spin-resolved 2p 3 / 2 
spectra, where the up- and down-spin electron removal 
excitations are considered separately. Because the SOI 
is not included in the calculation, we could not distin- 
guish 2p 3 / 2 and 2pi/ 2 . The dependence of the spectral 
shapes on the element, the excited core, and the spin 
of the core hole left behind are well reproduced by the 
calculation. The Fe 2p spectra show almost a single- 
peak structure for the up-spin and down-spin electron 
removals. In contrast, the spectral intensities for Ni are 
distributed in a wide range of binding energy with the 
notable satellite structure for the up-spin electron re- 
moval and the shoulder-like structure for the down-spin 
electron removal. The spectra for Co exhibit interme- 
diate features between those for Fe and Ni; the spec- 
tra are widely distributed around the higher binding en- 
ergy region with noticeable shoulder structure for the up- 
spin electron removal, while the spectra exhibit almost a 
single-peak structure for the down-spin electron removal. 
The satellite intensities in Ni 2p spectra are interpreted as 
coming from the final state that an electron on the qua- 
sibound state is excited to the unoccupied one-electron 
states on \f (0)). Such satellite or shoulder structures ex- 
ist in all elements Fe, Co, and Ni for the up-spin 3s elec- 
tron removal. 16 We have discussed the difference between 
the 2p and 3s spectra as well as the difference between 
Fe and Ni spectra in connection with the one-electron 
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Figure 5: Ni 2p XPS spectra. Gray thick curve represents 
the total intensity. The curves denoted / , up /' , dn / , 
/ , and _p 3 ' are the contribution from the final states in- 
cluding zero, one up-spin e-h pair, one down-spin e-h pair, 
two e-h pairs, and three e-h pairs on the final state |/(0)), 
respectively. Spin-resolved experimental data are taken from 
Ref. ED- 
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Figure 6: Difference of the threshold energy between the 
states that an up-spin core-electron is removed and that a 
down-spin core-electron is removed. Left, down, and right 
triangles represent Aw" 1 for the 2s, 2p, and 3s core-ionized 
states. Circle represents the local spin moment at the core 
hole site in /jb- Dashed lines are for a guide to eyes. 



states screening the core hole. We have explained the 
origin of these differences by a combined effect of the 
different magnitude of the 2p-Zd and 3s-3<i exchange po- 
tentials and the different occupation numbers in the 3d 
states in Fe and Ni. 

Although the calculation reproduce consistently the 
spectra shape depending on the elements, the excited 
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Figure 7: Difference of the up-spin and the down-spin core- 
level energies. Left, down, and right triangles represent Ae c 
for 2s, 2p, and 3s core-ionized states. Circle represents the 
local spin moment at the core hole site in fiB- Dashed lines 
are for a guide to eyes. 



core, and the spin of the core hole left behind, there is 
a clear discrepancy on the satellite position in Ni; 6eV 
satellite in Ni that the calculation gives the satellite is 
given at merely 4 eV higher than the threshold. One rea- 
son for this discrepancy may be that the fully screened 
potential approximation is not appropriate due to the 
insufficient time of screening core hole in Ni in compar- 
ison with the core hole lifetime. Another reason might 
be caused by the approximate nature of LDA. To clarify 
this issue, we need further studies. 

In connection with the magnetic splittings of the 
threshold energy, we have directly evaluated these val- 
ues within the fully screened potential approximation by 
using the KKR-CPA method. Olovsson et alM^HL have 
discussed the core-level chemical shifts in several metallic 
alloys using the same method, and pointed out that the 
fully screened potential approximation gives good cor- 
respondence to the experiments. The core-level energy 
depends on its spin, which difference is roughly propor- 
tional to the local spin moment. We have found that the 
magnetic splitting of the threshold energy decreases more 
rapidly than the local spin moment with moving from Fe 
to Ni; it is almost zero for Ni for the 2s, 2p, and 3s elec- 
tron removals, in spite of a finite local moment 0.6 /ib- 
The splitting may be directly measured by the shift of 
the peak around the threshold in the spin resolved 2s 
or 3s spectra. Actually, in the recent experiment of the 
2p 3 / 2 HAXPES spectra, 9 such a splitting is measured as 
0.9 eV. Our calculation also suggests that the splitting 
would hardly be observed in Ni. To be more precise, it 
may be better to take account of the SOI and to increase 
the super-cell size. 

We have considered only the one-electron states on the 
r point in the first Brillouin zone onto which electrons are 
distributed in the calculation of the XPS spectra. This 
may not cause large errors except the intensity near the 
threshold, because the first Brillouin zone is reduced to 
a smaller size in a system of larger supercells. With in- 
creasing the states onto which electrons are distributed, 
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we expect that the overlap between the lowest energy 
state in the presence of core hole and the ground states 
would be reduced, and that the contributions from e-h 
pair creation would increase near the threshold, leading 
to an asymmetric peak near the threshold as a function 
of binding energy. Such behavior has been demonstrated 
in numerical calculations on finite-size systems^i On 
the other hand, the structures in the high binding en- 
ergy region are expected to be only a little influenced by 
such a refined treatment. In any case, to be more quan- 
titative, we need to increase the size of supercells in the 
calculation. 

Another ab initio approach has been tried by using 
multiple scattering theory. 33 Since the experimental data 
have been accumulated for XPS spectra and the x-ray 



absorption spectra near the L edge, the extension of the 
present method to calculate the absorption spectra is left 
in future study. 
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